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The numerical determination of combustor flowflelds Is of great value to 
the combustor designer. An a priori knowledge of the flow behavior can speed 
the combustor design process and reduce the number of experimental test rigs 
required to arrive at an optimal design. Even 2-D steady Incompressible Iso- 
thermal flow predictions are of use; many codes of this kind are available 
[1], each employing different techniques to surmount the difficulties arising 
from the nonlinearity of the governing equations and from typically Irregular 
combustor geometries. We will look at mapping techniques (algebraic and ellip- 
tic PDE), and at adaptive grid methods (both multi-grid and grid embedding) as 
applied to axlsymmetrlc combustors. 
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ALGEBRAIC GRID GENERATION IN THE STARPIC CODE • 

t 

1 

The solution of the Turbulent Reynolds Equations for arbitrary geometries i 

can be handled In several ways [1]. One popular technique Is to represent a , 

boundary with discrete 90° steps. Atlernatlvely, the geometry can be mapped • 

Into a square solution plane via a variety of available transformations. The ’ 

simplest transformation Is an algebraic one as shown In Figure 2. i J 

On transforming the governing equations, new chain rule terms are Intro- , 

duced, each of which contains a dh/dx term which comes Into play only at axial 
locations where there Is a sloping or curved wall segment, thereby mathemati- 
cally affecting expansion and recirculation In the flow. The complexity of ! 

extra terms Is offset by the ability of the technique to handle general bound- 
ary shapes, as well as the simplification of boundary conditions In the mapped 
plane. For these reasons algebraic mapping shows great promise. Figures 1 | 

through 6 show the effects of algebraic mapping on the STARPIC code. • 

; j 

■ i 

ELLIPTIC PDE GRID GENERATION FOR A BIFURCATED j 

COMBUSTOR INLET DIFFUSER j 

The purpose of a combustor Inlet diffuser Is to convert kinetic energy to I 

pressure; a reliable prediction method for diffuser flows will lead to the 
efficacious design of those diffusers. In this work body fitted coordinate ! 

transformations are employed In the solution of turbulent flows In a bifur- 
cated combustor Inlet diffuser. The work Is aimed at comparing the numerical • : 

solution for different kinds of transformations with experimental results. 

Furthermore, different numerical methods will also be employed for comparison. ; 

Figures 7 through 12 Illustrate this work. 
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MULTIGRID AND GRID EMBEDDING TECHNIQUE FOR 
TURBULENT FLOWS IN COMBUSTOR 

A multigrid method [2] applied to axlsymmetrlc turbulent flows Is useful 
In speeding convergence of numerical schemes. This method reduces global 
errors by economical relaxation of errors on coarse meshes Instead of by labor 
Intensive relaxation on one fine mesh. The multigrid equations, a summary of 
the multigrid algorithm, and relationships between coarse and fine meshes for 
a staggered mesh are depicted In Figures 13 through 17. 

Grid embedding [3,4] similar to the multigrid method, uses a fine mesh In 
steep gradient regions to refine the solution locally, where needed. It elim- 
inates the wasteful process of using a fine grid globally. Including In regions 
where It Is not needed. Grid embedding Is depicted for a staggered mesh. As 
Indicated In the last figure, this method Is still under development. Like the 
multigrid method, It will steadily be extended to more arbitrary combustor 
geometries . 
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ADAPTED FROM 
STARPIC 


CONTRASTING STARP1C AND STARMAP 

ADVANTAGES OF MAPPING; 

ANY SINGLE VALUED FUNCTION VALID FOR h(,). 

SOLUTION OCCURS ON A SQUARE MESH IN THE < - n PLANE. 
BOUNDARY CONDITIONS ARE SIMPLIFIED. 

DISADVANTAGES; 

ADDITIONAL TERMS ARISE VIA TRANSFORMATION. 

HANDLING OF “WALL FUNCTIONS" NOT OBVIOUS. 
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SIAM 

NEARING OPERATION OF COMPLETELY MAPPED CODE. 

FOLLOWING FIRST SUCCESSFUL RUN, CODE WILL BE 
TESTED WITH A VARIETY OF GEOMETRIES IN TWO 
PHASES: 

PHASE (I) VALIDATION OF STARMAP IN 

DIRECT COMPARISON WITH STARPIC 
RESULTS; 

PHASE (II) GENERALIZATION TO OTHER GEOMETRIES 
NOT PRESENTLY HANDLED BY STARPIC. 
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ORIGINAL PAGE 13 
OF POOR QUALITY 
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STATUS 

TRANSFORMATION OF THE GEOMETRY HAS BEEN ACCOMPLISHED, 
TRANSFORMED GOVERNING EQUATIONS HAVE BEEN DISCRETIZED. 
COMPUTER CODES TO SOLVE TRANSFORMED EQUATIONS ARE 
BEING PREPARED. 
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MULTIGRID ALGORI THM 
GRID HIERARCHY 


G K " 2 

FUNDAMENTAL RELATIONS 

LET o*(r.x) - s <r.x> • SET OF TURBULENT REYNOLDS 

EQUATIONS 

IN OPERATOR NOTATION THIS BECOMES 
LU) * o 

FINITE DIFFERENCING ON GRID G K GIVES 

i k U k ) • ( k 

WHERE F K IS MADE IDENTICALLY ZERO ONLY ON THE 
FINEST MESH 
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CHECK 


ORIGINAL PAGE IE 
OF POOR QUALITY 
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IF THIS RATIO 
RELAX THE SET 

WHERE 


■0.9, SWITCH TO COARSER GRID. ON GRID G K_1 
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THUS. f k '’ CHANGES UPON DIFFERENT VISITS TO G*" 1 . 

ONCE * k '’ CONVERGES ON G K_1 , INTERPOLATION CAN SUPPLY 
MISSING VALUES IN G K . I.E.. 


♦ k ■ *Li 


WHERE ik., IS AN INTERPOLATION OPERATOR FOR 

g k * 1 to g k . 

FOR SMOOTHER RESULTS. A MODIFIED VERSION IS USED 
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IS THE INJECTION OPERATOR FROM GRID G K_1 
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HULTIGHID AND CRIP flUEDDIRG TECHH10IJES 
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V VELOCITIES ONLY, 
SHOWN FOR SIMPLICITY 


FLOW DOMAIN SUBDIVIDED INTO COARSE MESH REGION AND FINE MESH REGION. WHERE 
STEEP GRADIENTS REQUIRE MORE RESOLUTION; 

A’S ARE FINEST GRID FOR COARSE MESH REGION - IN FINE MESH REGION, FINE 
GRID IS USED TO REFINE SOLUTION AT A'S; 


BOUNDARY VALUES ARE UPDATED FOR EACH COARSE MESH SWEEP BY INTERPOLATION 
BETWEEN A 'S ON BOUNDARY; 

L k * k • F k in COARSE MESH REGION, AS BEFORE, AT A 'S; 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


STATUS. 

MULTIGRID METHOD IS BEING APPLIED TO THE STARPIC 
CODE 


FUTURE DEVELQfMS. 

REFINE EXISTING MULTIGRID METHOD 

DEVELOP EMBEDDED GRID METHOD 

EXTEND BOTH METHODS TO MORE ARBITRARY GEOMETRIES 
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